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1 Introduction 



There exist significant nuclear effects in the cases of J/^f, and even of Drell-Yan pair 
production. These effects are discussed in many papers (see below). We present the 
possibility to estimate them directly from the experimental data. 
Let us consider the variable 

P 

xo = —, (1) 
Po 

where p is the momentum of produced secondary and po is the initial momentum of the 
beam particle both in cm. frame. At high energies in the case of nucleon target the 
maximum value of x is close to unity ((x ) ma _ x — > 1). In the case of nuclear target the 
situation is more complicated because there are many different contributions. First of 
all there exists the so called cumulative production (with x > 1) PQI21E1- However it is 
a special process which we are not going to consider in this paper. 

Many other processes result in (xo) max < 1 due to the initial- and/or final-state in- 
teractions jlj. For example in the case of Drell-Yan production there exist initial-state 
effect such as energy loss of the incident quark in nuclear target as well as nuclear shad- 
owing In the case of J/^f production there is an additional source of suppression 
connected with final-state interactions 

Moreover, in all processes on nuclear targets (except of coherent ones) some fraction 
of energy is used for nuclear disintegration. The nuclear target is destroyed and several 
nucleons, as well as light nuclear fragments, (say, a-particles) appear in the final state. It 
can be considered as a phase space limitation. This fraction is numerically not so small 
(it is many times more than the nuclear binding energy) (HI Ej • This source of energy loss 
also leads to decrease of (xo) m3fX but it was not taken into account in all modern papers. 

In the present paper we will not consider every source of suppression separately, but 
we are going to consider them all together including the last one. 

It is clear that all initial-state energy losses are equivalent to decrease of the incident 
beam momentum p . As the phase space effects as well as final state interactions decrease 
(^o)max, we assume that their influence in the considered processes can be effectively 
described by the same way, as the additional decrease of p . So we change p by p Q — pa, 
where pa is a phenomenological parameter which accounts for all energy losses. 

In what follows we will consider the ^-dependences of J/ty and Drell-Yan pair pro- 
duction in terms of Xf: 

x F = (2) 

PO -PA 

and we assume that it is possible to find the shift pa from the condition that the ratio 
of multiplicities on nucleon and nuclear targets 

RhA/h P (x F ) = const(x F ) ~ 1, (3) 
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whereas the same multiplicity ratio in terms of Xq has evident A-dependence 



RhA/h P (x ) = f(x ). 



(4) 



Evidently, such rescale is reasonable only for not very small x$ values. 

We will determine the shift pa from the experimental data and we will compare 
them with several independent estimations. Such approach allows us, in particular, to 
investigate the energy dependence of all nuclear effects. In conclusion we will compare 
our results with theoretical calculations E] • 

2 A-dependence of J /iff production at large 

Charmonium production off nuclei has drawn much attention during the last two decades, 
since the NA3 experiment at CERN jTUJ has found a steep increase of nuclear suppression 
with rising Xq. This effect was confirmed later in the same energy range |TT] as well as 
at high energies [121 HSl El • 

For the purpose of our analysis we used the experimental data on J '/ty production in 
proton- nucleus (pA) collisions published in [T21 EH US] • We also used data on J/ty pro- 
duction in n~A collisions [TTJ [THJ - We analysed the ratio Ra 1 /a 2 °f inclusive differential 
cross sections for J '/\& hadroproduction on A\ nucleus to that on A 2 : 



The analysis was perfomed in the following way. Suppose one has two xo-spectra for 
jj-fy production on nucleon and nuclear targets or on two different nuclei (light and heavy 
ones). Then one can shift the spectrum that corresponds to the heavy nucleus according 
to Eq. (J2J) by changing pa parameter. Assuming that nuclear effects are small for very 
light nuclei, it is possible then to calculate the fraction of beam energy/momentum spent 
on nuclear effects. This may be done by calculating the ratio of the shifted spectrum to 
the spectrum that corresponds to the light nucleus. When this ratio is close to unity then 
the corresponding shift will give the absolute value of energy /momentum loss caused by 
the mentioned nuclear effects. 

Let us analyse the spectra [H] presented in Fig. El The spectra represent the differ- 
ential cross sections for J/^'s produced inclusively in 800 GeV/ c pCu and pBe collisions 
measured by E789 Collaboration. One can see the evident nuclear suppression for copper 
target. 

Now to calculate how much energy/momentum of the beam particle is spent on nuclear 
effects we should shift the spectrum for the heavy nucleus (namely copper). The shifted 
spectrum is presented on Fig. El The solid curve in Fig. |21 represents the fit to the shifted 
spectrum needed for further calculations of the ratio of the presented spectra. 




(5) 
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The calculated ratios Eq. (jSJ) one can see in Fig. EI Solid circles correspond to the 
original ratio measured by E789 Collaboration, two other ratios were calculated for 
different shifts: pa = 0.5 GeV/c (solid squares) and pa = 0.7 GeV/c (open squares). 
These pa values represent the amount of absolute momentum spent on nuclear effects. 

The values of pa in cm. frame correspond to the shift 

Ax c F m - = (6) 
Po 

For large xp Axp b ' ~ Axp m ', so we can calculate the absolute value of energy losses in 
lab. frame. In our case the shift pa = 0.5 GeV/c corresponds in the lab. frame to the 
energy losses Ap^u/Bc ~ 20 GeV/ c. 

Unfortunately E789 is the only experiment that measured absolute cross sections for 
J/ty production in pA collisions for both light and heavy nuclei. Experimentally it is 
much easier to measure ratios of the cross sections at once, thus most collaborations 
present only the ratios (without absolute values of the cross sections) as their results. 
Consequently it is difficult to analyse those experiments in such a way we did above. 

However the most recent and most precise experiment on J/ty production by E866 
Collaboration |T2j was performed at the same energy (800 GeV) as E789 and used Be 
and W targets. Thus using the E789 Xo-spectrum for beryllium target and E866 ratio of 
spectrum on beryllium to that on tungsten one can extract the absolute J/^> production 
cross section for tungsten. Since the Xo-scale covered by E866 experiment (—0.1 < 
x < 0.93) is larger than one covered by E789 (0.3 < x < 0.95), then to extract E866 
spectrum for W target we combined the E789 data with the data obtained by E672 and 
E706 Collaborations for Be target in the range of 0.0 < xq < 0.5 |15j . Fig.|3]represents the 
ratio measured by E866 (Fig. and the absolute spectra for Be target (solid squares) 
combined from the mentioned data sets, and for W target (open squares) extracted from 
the ratio (Fig. HJ>). Now it is possible to analyse the last two spectra in the same way as 
was done in the case of E789 data (see Fig. EJ). We omit intermediate calculations and 
present the final result. The shifted ratios at two different shift values (pa = 1.2 GeV/c 
- solid squares, and pa = 1-5 GeV/ c — open squares) are presented in Fig. El This 
corresponds to the absolute energy losses in the lab. frame Ap l ^ Be « 50 GeV/ c. 

One can see that the absolute value of energy /momentum loss in tungsten is more 
than two times larger than one in copper target, which is rather clear. Indeed, W target 
is ~ 2.9 times heavier than Cu. Consequently, nuclear effects in the former should be 
stronger, thus the value of energy/momentum lost by the projectile in W target should 
be larger than that in Cu. 

Besides proton induced reactions we considered the data on J/ty production in tt~A 
collisions. The data available are the ratio of differential cross sections for J/\I/ production 
on W to that on Be target at 125 GeV/c (TTj (Fig. Ek)> and absolute differential cross 
section for J/ty production on Be target at 150 GeV/c [TH]- Since the energies for the 
two data sets slightly differ from each other, we build the corresponding invariant cross 
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section from that in Ref. ^U], which is presented in Fig. Eb- To extract the spectrum for 
a heavy nucleus (namely, tungsten) we applied the procedure described above. Omitting 
intermediate results we present the range of pa values obtained from the analysis (Fig.|7|): 

p A = 1.5 - 1.7 GeV/c, (7) 

which correspond to Ap l ^ Bc m 12 GeV/c. The obtained Ap l ^ Be value for ir~A collisions 
at 125 GeV/ c is about 4 times smaller than that obtained for pA collisions at 800 GeV/ c. 
Some part of this difference (say, factor ~ 1.5) can be connected with smaller pion- 
nucleon (irN) cross section in comparison with NN cross section. Another part of the 
difference can be connected with the dependence of nuclear effects on the initial energy. 

Unfortunately the errors of the data combined with the analysis errors result in too 
large error of the final result, which does not allow us do draw a definite conclusion. 



3 ^-dependence of Drell-Yan production at large xp 

Since the Drell-Yan mechanism produces lepton pairs which only interact electromagnet- 
ically, the ^4-dependence is expected to be weak because no final-state interactions affect 
the lepton pair. However some initial-state interactions may affect the A dependence. 

Almost all the experimental results on Drell-Yan production are presented in terms 
of the ratio of inclusive differential cross sections on a heavy nucleus to that on a light 
one. And there is no opportunity to extract desired spectra separately as was done in 
the previous section. 

However we developed a Monte-Carlo (MC) event generator (HARDPING — Hard 
Probe Interaction Generator) that extends well known HIJING MC jTTj on Drell-Yan 
pair production process and some initial-state effects are accounted for [THJ EH] . HARD- 
PING MC describes well the data on Drell-Yan pair production in hadron-nucleus (hA) 
collisions at high-energies [TH] ITU] . 

Using HARDPING MC we simulated absolute spectra for Drell-Yan pair production 
in pA collisions on W and Be targets at 800 GeV/ c. Then we applied the procedure 
described above for the two simulated spectra. The results of the analysis are presented 
in Fig. |H| The figure represents the original E866 data jS] (solid circles), simulated ratio 
without any shift (open circles) to demonstrate consistency between the simulated results 
and the experimental data, and shifted ratios obtained with HARDPING MC. As was 
predicted the fraction of energy /momentum lost by the projectile on nuclear effects is 
small for the case of Drell-Yan pair production. This is because there is no final-state 
interactions in Drell-Yan production process. 

The obtained value for pa ~ 0.3 GeV/c corresponds to Ap^ b ; Bo ps 12 GeV/c. 
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4 Conclusion 



In summary, we considered the energy/momentum losses of the projectile in hA collisions 
at high-energies from the available experimental data. What we were interested in is how 
much energy /momentum of the projectile is spent on all the nuclear effects including the 
effect of nuclear disintegration. 

The energy losses estimated from the experimental data are in reasonable agreement 
with theoretical calculations [HI Ej. Namely, our result Ap^ b y Bc pa 50 GeV/c for J/\l/ 
production at 800 GeV/c correspond to energy loss rate 

dE/dz w 5 GeV/fm, (8) 

(We assume the length of full trajectory ~ 1.5Ra and we neglect nuclear effects in Be 
target). The analysis 7| predicts dE/dz ~ 3 GeV/fm for initial-state quark energy 
loss rate, i.e. energy loss before the hard interaction point. However initial-state quark 
energy losses is not the dominant effect in J/ty production processes on nuclear targets. 
The main contribution to the suppression of J/^'s arises from final-state interactions. 
Also there exists strong gluon shadowing at large Xq as well as gluon enhancement at 
small xq. Nevertheless, Ref. [7j does not present how much projectile energy is spent on 
each nuclear effect. (The effects of gluon shadowing for J/ty production were calculated 
in Ref. |20j.) Thus we assume that main part of our result Eq. (JSJ) is explained by the 
mentioned effects. 

The estimations of energy losses in pairs Be-W and Be-Cu are larger than the ratio 
of A 1 / 3 values (i.e. length of trajectory) for W and Cu nuclei. It can be connected with 
rather large error bars, or with the A-dependence of energy losses more strong than A 1 / 3 
behavior. The last reason is not excluded because the multiplicity of secondary protons 
produced in hadron-nucleus collisions in the nuclear targets fragmentation region with 
energies < 1 GeV has A-dependence more strong than A 1 ^ 3 \'21\ 12*2] . On the other hand, 
the energies of these protons are determined by energy losses of the incident particle. 

In the case of Drell-Yan pair production calculations jHj predict quark energy loss rate 
dE/dz pa 3 GeV/fm. There also exists shadowing effect in nuclear target, however there 
is no numerical estimate for this effect in Ref. jHj, it was considered in Ref. "23|. Our 
result Ap l ^ Be pa 12 GeV/c obtained for Drell-Yan production corresponds to energy 
loss rate dE/dz pa 1.2 GeV/fm for full trajectory. However one should take into account 
the path of the projectile before the hard interaction point only, which is ~ 4.4 fm for 
W target 6 . Thus dE/dz pa 2.7 GeV/fm. Assuming that shadowing effect is small 
for Drell-Yan production at 800 GeV/ c we conlude that our results are in reasonable 
agreement with the ["J] calculations. 

It is necessary to note that some part of QCD energy losses can be used for nuclear 
destruction and fragmentation due to the final-state interactions. This can explain rather 
large energy of secondary target nucleons and nuclear fragments observed in [HI EI - 

From the difference in the values of Ap^!y Be for pA and it A collisions we can conclude 
that nuclear effects probably depend on projectile energy. 
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FNAL E789 

o p + Cu -> J/\|/ at 800 GeV/c 
• p + Be -» J/\|/ at 800 GeV/c 
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Figure 1: Differential cross section for J/ty production in pCu and pBe collisions at 
800 GeV/c 
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o p + Cu -» J/\j/ (shifted) 
• p + Be — > J/\|/ 
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Figure 2: The same as in Fig.^ but the spectrum for Cu was shifted according to Eq. 
with p A = 0.5 GeV/c. 
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Figure 3: The ratios of inclusive differential cross sections, calculated for different pa 
values. 
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(a) 

pA -> J/\|/ at 800 GeV/c (E866/NuSea) 
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Figure 4: (a) The ratio of inclusive differential cross section for J/\& production on W 
target to that on Be measured by E866 [T^]. (b) Combined spectrum (see text) for Be 
target (solid squares), and spectrum for W (open squares) target extracted from the 
ratio shown in Fig. 
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Figure 5: The ratios of inclusive differential cross sections, calculated for different pa 
values. 
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(b) 
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Invariant cross-section for: 
TV + Be -> J/\|/ at 150 GeV/c 
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Figure 6: xo-spectra for J/^'s produced in ir~A collisions, (a) The ratios of inclusive 
differential cross sections for J/ty production on tungsten to that on beryllium [11J. (b) 
The invariant differential cross section for J/ty production on Be target. 
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Figure 7: The ratios of inclusive differential cross sections for J/ty production in n A 
collisions at 125 GeV/ c, calculated for different pa values. 
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Figure 8: The ratios of inclusive differential cross sections for Drell-Yan pair production 
in pA collisions at 800 GeV/c, calculated for different pa values. 
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